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UNITED KINGDOM 
ABSTRACT 

Recently there has been much work devoted to considering some of the many and varied 
interaction mechanisms which may be operative in three-dimensional boundary layer flows. 
Here we are concerned with resonant triads of crossflow vortices. In contrast to much of the 
previous work we examine the effects of interactions upon resonant triads where each member 
of the triad has the property of being linearly neutrally stable; then the importance of the 
interplay between modes can be relatively easily assessed. We concentrate on investigating 
modes within the boundary layer flow above a rotating disc; this choice is motivated by 
the similarity between this disc flow and many important practical flows and, secondly, our 
selected flow is am exact solution of the Navier-Stokes equations which makes its theoretical 
analysis especially attractive. Firstly we demonstrate that the desired triads of linearly 
neutrally stable modes can exist within the chosen boundary layer flow and then subsequently 
obtain evolution equations to describe the development of the amplitudes of these modes once 
the interaction mechanism is accounted for. It is found that the coefficients of the interaction 
terms within the evolution equations are, in general, given by quite intricate expressions 
although some elementary numerical work shows that the evaluation of these coefficients 
is practicable. The basis of our work lends itself to generalisation to more complicated 
boundary layers and effects of detuning or non-parallelism could be provided for within the 
asymptotic framework. 


1 Research was supported by the National Aeronautics and Space Administration under NASA Contract 
No. NAS1-18605 while the authors were in residence at the Institute for Computer Applications in Science 
and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 
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§1. Introduction. 


The desire to understand the phenomenon of transition to turbulence in fluid flows 
has motivated a wealth of studies concerned with the instability of three-dimensional 
boundary layer flows. Here we address the problem of providing a self-consistent 
asymptotic description of the interactions within a resonant triad of travelling cross- 
flow vortices occurring in a three-dimensional boundary layer at a large Reynolds 
number. In particular, we relate our analysis to vortices in the boundary layer of the 
flow induced by a rotating disc. This flow is susceptible to instabilities similar to those 
which occur in flows over a swept wing; a situation which has practical relevance to 
the development of Laminar Flow Control wings. Further, there is an exact solution 
of the Navier Stokes equations which describes the rotating disc flow which makes 
this particular flow especially suitable for a theoretical analysis of the interaction of 
travelling cross-flow vortices. 

The resonant interaction of a triad of waves is one of a number of possible mecha- 
nisms that play important roles in describing the nonlinear evolution of disturbances. 
Such an interacting triple is possible whenever three waves with wavenumbers k, and 
frequencies f lj (j = 1,2,3) such that ki = k 2 + k 3 and flj = fl2 + ^3 can coexist. 
This type of mechanism is the subject of our current investigation, but an excellent 
account of the variety of interaction phenomena which are possible within the scenario 
of fluid mechanics has been given by Craik (1). 

The crossflow vortex instability structure considered here can only occur in three- 
dimensional boundary layers and was first examined both theoretically and experimen- 
tally by Gregory, Stuart & Walker (2) in the context of rotating disc flow. They showed 
with the china-clay technique that this instability takes the form of a regularly spaced 
pattern of equiangular spiral vortices which are stationary relative to the disc. In this 
paper GSW showed, using inviscid theory, that the stationary vortex mechanism is as- 
sociated with an inflexional profile in which the inflexion point coincides with a point 
of zero velocity somewhere in the flow. Their calculation of the number of vortices 
expected to be seen in the flow above a rotating disc gave a number approximately 
four times greater than the observed value of about 30, but the angle of 13° between 
the axes of the vortices and the radius vector predicted was in good agreement with 
that seen in experiments. Brown (3), and later Cebeci fc Stewartson (4) solved the 
Orr-Sommerfeld equation numerically and obtained values for the critical Reynolds 
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number for transition much lower than that observed by Gregory et. al. and used 
by Stuart in his calculation for the number of vortices. Malik (5) found the neu- 
tral stability curve for these stationary disturbances and he showed the existence of 
a second stationary mode (this one however being viscous) which corresponds to zero 
wall shear stress of the crossflow velocity profile. The asymptotic structures (at high 
Reynolds numbers) of both this viscous mode and of the inviscid mode of GSW were 
obtained by Hall (6). The linear analysis presented in that paper has been extended 
by MacKerrell (7) and Gajjar (8) to include nonlinear effects. 

The aforementioned investigations have concentrated on stationary modes. How- 
ever, Faller & Kaylor (9) and Lilly (10), in their studies of periodic perturbations 
imposed on Ekman boundary layer flows, mentioned the existence of an instability 
with non— zero phase speed. Federov et. al. (11) showed that in rotating disc flow 
a similar instability occurs much earlier than the stationary mode observed in the 
experiments of Gregory et. al. This instability also appeared as a pattern of spiral 
vortices, but which was moving relative to the disc. The number of these vortices was 
observed to be between 14 and 16, and had axes inclined at about 20° to the radius 
vector. Recently, Bassom & Gajjar (12) have conducted a theoretical analysis of some 
aspects of non-stationary vortices in rotating disc flow and, in particular, investigated 
the properties of neutrally stable nonlinear modes. Our present concern is with the 
time dependent version of the high Reynolds number inviscid crossflow vortex whose 
flow structure was elucidated by Hall (6). The modifications required to his anal- 
ysis to incorporate the effects of unsteadiness are relatively straightforward and are 
considered in the following section. 

There has been much work in relation to the important problem of the interaction 
of Tollmien-Schlichting disturbances with Gortler vortices (an instability associated 
with flows over curved surfaces). See, for example, (13-17) and the references therein. 
Also, some study has been made of the interaction between Tollmien-Schlichting waves 
and steady crossflow vortices, (18-20), and of the triad interaction between three un- 
steady crossflow vortices under consideration here. In particular, we refer to the work 
of Reed (18, 19) and El-Hady (21), who found that the growth of each of the insta- 
bility modes could be substantially influenced by their mutual interaction mechanisms. 
Additionally, non-parallel effects were included and it was concluded that the evolu- 
tion of the vortices was critically dependent upon the initial amplitudes of the modes. 
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However, these accounts ignored the crucial fact that for a completely rational descrip- 
tion of the importance of the interaction process in the development of the disturbances 
each of the instability modes involved should be neutrally stable at leading order in 
their own rights. Otherwise, the vortex growth due to the interaction mechanism is 
no larger than would be experienced by the modes in the absence of the interaction, 
and the true importance of the mechanism is extremely difficult, if not impossible, to 
accurately assess. For this reason, our present aim is to first illustrate that a suitable 
triad of travelling crossflow vortices can exist in a three-dimensional boundary layer; 
crucially, these modes are chosen so as to ensure that each is individually neutrally 
stable at leading order. 

Once the existence of the desired triad is established we move on to consider the 
problem when all three modes are simultaneously present. In the case of a single, 
infinitesimally small mode the asymptotic flow structure consists of a main layer (co- 
incident with the majority of the boundary layer) augmented by a thin, viscous wall 
layer next to the surface of the disc and by at least one linear critical layer present 
somewhere in the flow. With all modes of the triad present, the corresponding flow 
structure becomes much more complicated, for now we have at least three critical lay- 
ers (one for each mode) and sometimes more. It transpires that in order to obtain a 
rational description of the evolution of the triad each of the critical layers needs to be 
carefully analysed, for each contributes non-zero terms to each of the final, coupled 
amplitude equations for the three modes. 

An important characteristic of our analysis is that each vortex is taken to be suffi- 
ciently small and to develop on appropriately slow length and timescales such that all 
the critical layers remain linear in character. We emphasize that although the equa- 
tions obtained in this paper are specifically oriented towards the rotating disc flow, for 
the cases of other basic velocity profiles the derivation procedure would be essentially 
identical. The method adopted also allows nonparallel effects to be accounted for with- 
out formal difficulty. The linked triad of evolution equations, somewhat unfortunately, 
suffer from the handicap that they are very complicated and require much delicate non- 
trivial numerical work to establish the interaction coefficients for any particular flow. 
Here we perform some of the calculations for the rotating disc situation to demonstrate 
that the derivation of the evolution equations is practicable. The equations so derived 
have complex-valued coefficients and are not dissimilar in form to those derived by 
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Smith & Stewart (22) who were concerned with high frequency resonant triads within 
boundary layers. We consider that the principal result of the current work is that in 
three-dimensional boundary layers resonant triads of non-stationary crossflow vortices 
do exist. Further, given such a triad of modes which are neutrally stable at leading or- 
der, we show that we can develop a fully rational asymptotic account of the interaction 
processes. 

The procedure for the remainder of the paper is as follows. In section 2 we 
briefly consider the asymptotic structure of a linear, travelling crossflow vortex in a 
three-dimensional boundary layer, making specific reference to the case of rotating 
disc flow as an example. Numerical work in section 3 establishes that the desired form 
of a resonant triad of modes exists for this special basic flow and, armed with this 
information, we make a careful study of the interaction mechanism in section 4 to 
enable us to derive a coupled triple of evolution equations for the modes. Finally we 
perform some numerical work to establish the interaction coefficients, discuss a few 
solutions of the relevant evolution equations and we conclude with some discussion 
and suggestions for further work. 


§2. The basic equations and flow structure for a linear crossflow vortex. 

We have already stated that we intend to make an investigation of the interaction 
processes which take place between the three members of a resonant triad of crossflow 
vort ices and a natural first step is to consider the distu^ance structure for each indi- 
vidual mode in isolation. Here we concentrate on the rotating disc flow problem for 
definiteness, although the work can be trivially modified to account for other three- 
dimensional basic flows; for example the asymptotic suction profile or the attachment 
line flow, see Hall & MacKerrell (23). In the particular problem at hand, we consider 
the case in which the disc rotates about the 2 -axis with angular velocity 0. Relative 
to cylindrical polar co-ordinates (r, 6 , z) which rotate with the disc and in which r and 
2 have been made dimensionless with some reference lengthscale L , the continuity and 
Navier-Stokes equations for an incompressible fluid in the region z > 0 are 

V.u = 0, 

du 1 

— + (u.V) u 2 (ft A u) + ft A (ft A r) = Vp + z/V 2 u. 

at p 
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(2.1a) 

(2.16) 


Here u is the velocity vector, r is the co-ordinate vector, p is the fluid density, p is the 
fluid pressure and u the kinematic viscosity. The Reynolds number R for the flow is 
given by R = HL 2 /v and is taken to be large throughout the following. 

It is found convenient to define the small parameter e = R~ • . As the axes rotate 
with the disc, the basic flow is given by the Von-Kdrmdn solution 

u = u B = Ln(rfl(0,rP(0,e 3 ®(0) » p = pL 2 Sl 2 e 6 p( 0, (2.2) 

where 2 = e 3 £ and 0,0, © satisfy the equations 

0 2 -(l + V ) 2 + u‘xD = n", 20(1 + 0) + 0'© = 0", (2.3a, 6) 

©'+2G = 0. (2.3c) 

Here primes denote differentiations with respect to £ and the appropriate boundary 
conditions are 

0 = 0 = © = O on £ = 0; 0 — > 0, 0 — > —1 as £ — > 00 . 

The structure of an infinitesimal unsteady crossflow vortex imposed on the above 
basic flow follows closely that given by Hall (6) for the steady case. The mode has 
wavelengths scaled on the boundary layer thickness and an 0(1) wavespeed and so we 
consider perturbation quantities with r and 0 dependence given by E, defined by 

E = exp -3 (/ ot(r,e)dr + 00(e) - tfl(e)) j . (2.4a) 

We shall restrict our attention to examining disturbances in a neighbourhood of 
some point r = r„ and expand the wavenumbers and frequency in the forms 

a — ao + eai + . . . , /? = Po + e/?l + • • • > fl = Ho + + (2.46) 

In general these quantities will be complex but here we shall be concentrating on 
neutrally stable disturbances (at least to leading order) and so we take ao, /?o an <l ^0 
to be real. The disturbance structure in the z — direction is sketched in figure (1). There 
is an inviscid layer (I) of thickness 0(e 3 ) (the same depth as the boundary layer), and 
to satisfy no-slip conditions at the wall a viscous layer (II) of thickness O (e^ 2 ) must 
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be present. If the basic flow is perturbed so that the total non-dimensional flow field 
is given by 

( u , v, w,p ) = (rU,rV, e 3 tz>, e 6 p) + 6 [(17, V, W\P) E + c.c .] + . . . , (2.5a) 

where { < 1, then in zone (I) the disturbance quantities U, V, W, P (which are func- 

tions of £ alone) expand in the forms 

(U, V,W, P) = (Uo + eUl + ..., V 0 + eVi + ..., W 0 + eWi + ...,P 0 + ePi + ...). 

(2.56) 

Substituting expressions (2.5) in the continuity and the linearised Navier-Stokes 
equations yields at leading orders 

i ( a 0 U 0 + ^ Vo) +Wq = 0, (2.6a) 

' r„ / 

i (aor„Q + @oV — flo) 17o + r n U Wo = — iao-Po, (2.66) 

i (aor n U + PoV - Qo) Vo + r n V Wo = — — Po, (2.6c) 

and 

i (a 0 r n U + PqV — fl 0 ) JV 0 = -Po- (2.6d) 

If we define 

Uj = Qjr n U + PjV- 

then we obtain the classical Rayleigh equation 

tfo (w'o - 70 ^Vo) = K Wo, (2.7 a) 


where 7$ = • The required boundary conditions on (2.7a) are, firstly, 

|IV 0 (oo)| < 00, IV 0 (0+) = 0, (2.76, c) 

to avoid exponential growth at the edge of the boundary layer and to give the typical 
inviscid condition of tangential flow as the wall is approached. Technically, to reduce 
this slip velocity to zero at the surface of the disc, analysis of the viscous layer (II) is 
required but such analysis does not affect the leading order quantities c*o, Po and flo 
which are of the primary interest here. 
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The Rayleigh equation (2.7a) additionally has singularities at any positions where 
Uq = 0, and, if this occurs, the presence of an internal critical layer (at that position) 
is suggested. The details of the asymptotic analysis of such a (linear) critical layer 
follow well-known lines (see for example Drazin & Reid (24)) and so we merely sketch 
a few of the relevant results here. If in the vicinity of the critical layer, at £ = f say, 
we have the Taylor expansion 

o, = £%-(< -O'. ( 2 - 8a > 

where we have defined By* = (ayr n o(*) + /3y ?>(*)) |{=f* then the solution of (2.7a) as 
( — > £+ is given by 

Wo = X) 9r U ~ 0 r + E Gr U - 0 T In U ~ 0 , (2.8b) 

r=0 r=l 

where go,gi are constants upon which the remaining coefficients g n ( n = 2,3,...), 
G n ( n = 1,2,...) depend. In particular, we have G\ = <70-^02/ Boi- The solution 
(2.86) implies that Uq and Vq in (2.5) have logarithmic singularities as £ — ♦ £+ and 
to smooth these we need to invoke viscous effects in a thin layer ( III in figure 1) 
surrounding £ = £. Balancing inertial and diffusive terms suggests that this zone is 
of thickness 0(e 4 ) and if we define ( = £ + eip here, where ip = 0(1), the disturbance 


quantities develop as 

U = e~ 1 ^ 0 + In € &l + U 2 + eln + c ^4 + • ■ ■ > (2.9a) 

V = e ~ 1 Vo + In + f In CV 3 + tV* + • ■ ■ > (2.96) 

W = Wo + elnet^i + eW 2 + e 2 In eW 3 + e 2 W 4 + ..., (2.9c) 

P = Pq + e In ePi + eP 2 + e 2 In ePz + e 2 P4 + (2.9d) 


Briefly, we find on substituting (2.9) into the linearised Navier Stokes equations 

that 
V 0 = 

ao - Tj $r) Po { exp [- iTs9n(Btii) (* + f^) 1 exp ( - 3jkii) dT ' 

(2.10a) 
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where Po = Po(£)- This solution, as \tp\ — > oo, matches with the solution for Uq in 
(2.6) as £ — > £. Across the critical layer we also find that 


and 

d 2 W 2 iB 02 g 0 [°° 

diP 2 | -Soil K ^ 


Wi = 


w o = - 

g o-5q2 

B 01 


HnPo 

Boi ’ 

(2.106) 

( . B,n\ 


(2.10c) 

(. Bio'll { 

r 3 


3 |B 0 lK {2 ' m 


If we integrate (2.10d) once and then apply the asymptotic results given by Haber- 
man (25) concerning the behaviour of integrals of the general type 



dr, 


for large argument rp, we may show that the derivative dW 2 /dxp suffers a discontinuity 
across the critical layer and in particular that 


“ a T V' * OO 

dw 2 _ ITT Pp2g0 

- dip J ~ iPoil 

► - OO 


( 2 . 11 ) 


The implication of the solutions (2.10), (2.11) is that the Frobenius solution (2.86) 
for the eigenfunction Wq as £ — ► £+ continues to hold in the limit £ — > so 
long as the replacement of ln(£ - $) by (ln(£ - 0 - i 7 rsgn(B 01 )) in (2.86) is made. 
Equivalently, we have that as f , 


oo oo 

^0 = Eft((-(T+EG r ((-(T Me - 0 - iirsgn(B 01 )} , (2.12) 

r =0 r =1 

just below the critical layer. This is the well known phase shift of ‘iV’ across a linear 
critical layer. 

We now have a complete specification of the problem to determine the leading 
order neutral wavenumbers ao> @0 at some given frequency flo. The Rayleigh equation 
(2.7a) needs to be solved subject to the boundary conditions (2.76, c) and where the 
replacement given by (2.12) in the local series solution (2.86) across the critical layers 
is applied. The analysis presented here is specific to the case of rotating disc flow, 
but may easily be modified to account for a variety of other three-dimensional basic 
velocity profiles. 
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§3. The isolation of resonant triads. 


At the outset our aim was to verify that resonant triads (with the property that 
each member of the triad is neutrally stable at leading order) are possible within 
particular three-dimensional boundary layer flows. In the context of rotating disc 
flow, we therefore wish to seek ordered triples > ftoj >^0; ) 5 j = 1)2,3} such that 

(X oi = 002 + a 03> Pol = /?02 + /^03) ^01 = ^02 + ^03> (3-1) 

and with the characteristic that each triple is an eigenvalue set of the system (2.7), 
(2.8) k (2.12). To this end a numerical approach was adopted. 

A first concern is to find the location(s) of the critical layer(s). These occur 
wherever Vq = 0; i.e. wherever 

aor n U + (3qV = HO) (3.2a) 

or, if we define the leading order wavespeed co = Oo/r n ao; wherever 

U t = U + fxV = co, (3.26) 

with n = (3o/(r n ao). Hence, the ‘effective’ basic velocity profile U , is a function of fx 
and this flow is sketched in figure (2). It is known that a'(0) = 0.510, 0 (0) = -0.616 
and then as the parameter fx varies the effective basic velocity acquires a number of 
distinct possible forms. The boundary conditions in (2.3) suggest that U„ — * -fx as 

£ > oo and for /z > fx e (« 0.8284) the profile is a monotonic decreasing function. 

For 0 < (i < fi e , U, vanishes somewhere other than at the origin, and indeed when 
H « 0.235, the point of inflexion of the basic flow U, coincides with the point at 
which it vanishes. In this case the boundary layer becomes susceptible to the inviscid 
stationary vortex found by GSW. For n < 0, U t is always a positive function and for 
fj, —0.6 it is monotonic as well. We see from figure (2) that, depending upon the 
values of fi and co in (3.26), we may have either one or two critical layers. 

The numerical scheme used in practice requires that a value of fi be preselected 
and then ao> Po &n-d Oo are determined such that fx = Po/{r n ao). It is easily shown 
from (2.7a) that as f — > oo the eigenfunction Wo is proportional to exp(— 7of)- 
A suitably large value of £, say , was chosen and the asymptotic behaviour of 
the eigensolution was assumed to hold there. With these starting conditions, the 
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Rayleigh equation (2.7a) was integrated with a fourth order Runge-Kutta scheme to 
a point a small distance, say A, above the first critical layer encountered. There the 
series solution (2.86) was assumed to hold and go and g\ were chosen so that the 
necessary matching conditions held. Then, using the replacement suggested by (2.12), 
the solution just below the critical layer was established which could then be used to 
integrate (2.7a) until either the next critical layer position was reached (in which case 
the process described above was repeated) or until f = 0. In general, the predicted 
value of Wo(£ = 0) is complex and non-zero and so a Newton iteration on the real 
values ao and co was employed to ensure that the required condition (2.7c) was indeed 
met. Given the eigenvalues c*o and co for some chosen value jx it is then a simple 
task to retrieve the leading order azimuthal wavenumber flo and the frequency Oo- 
This numerical work was checked for independence from the chosen value , the step 
length used in the Runge-Kutta integration and of the distance A from the critical 
layer (s) at which the series solutions were applied. 

The solutions of the eigenproblem for various values of the parameter /x (which 
is tan" 1 (f), where 4> is the angle which the considered disturbances make with the 
outward radial direction on the disc) are shown in figure (3), where we illustrate the 
dependence of the leading order radial wavenumber ao and the frequency fto on the 
wavenumber (3q. There are several noteworthy features of the solutions. Firstly, there 
are two distinct branches in solution space which are independent of each other and do 
not coalesce. On one branch (labelled I on figure (3)) we see that there is a non-trivial 
solution of the eigensystem with zero frequency and this corresponds to the inviscid 
stationary mode of (2). Here ao 1 , 1 ? Po ?v 0i26r^ and jjl w 0.235^ and tli6ii the basic 
profile U t has a point of inflexion coincident with the point at which it vanishes. As 
the parameter /x increases we find that the frequency becomes negative and using the 
sketches of the basic flow in figure (2) we see that we are reduced to an eigenproblem 
with a single critical point. Appealing to the ideas of Stuart in (2), in the case of 
only one critical layer the neutral wavespeed co is determined by the condition that 
the function U, in (3.26) must have a point of inflexion where U t = co. 

For a waveangle parameter g. less than the inviscid stationary mode value of 0.235 
we now have two critical layers in the flow and the neutral frequency is positive. As 
H — * 0 the wavenumbers and frequency become small. Finally for this branch, as 
H — * 0 or n — > oo then the wavenumbers in the radial and azimuthal directions 
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become of different orders of magnitude and the symmetry of the governing equations 
(2.6) is lost. The structures which would govern these limiting cases are then not 
covered by the work of section 2 and it would be interesting to identify and analyse 
these structures. 


On the second branch of solutions (labelled II on figure 3) it is immediately seen 
that only positive frequencies are found. In particular, we remark that as /z increases 
and approaches ji c « 0.187 (corresponding to a disturbance waveangle of 10.6°) the 
wavenumbers and frequency become small. This suggests the existence of a longer 
wavelength disturbance with a small wavespeed and this structure corresponds to that 
of the ‘upper-branch’ mode investigated in (12). In this paper, small perturbations 
with wavenumbers O(R^) and wavespeeds 0(R~ ”) were shown to be possible at 
waveangles between approximately 10.6 s and 39.6 s . As /z decreases from /z e in our 
current problem and becomes negative the convergence of the numerical scheme em- 
ployed becomes increasingly difficult to achieve and only the neutral solutions which 
can be reliably located by our method have been shown in figure (3). Indeed, for 
sufficiently negative /z the basic profile is monotonically increasing and has no point 
of inflexion. At these values of parameter, no neutral solutions of the sought type can 
exist and so it suggests the existence of some /z„ < 0 at which our solutions cut out. It 
is seen from figure (2) that as n decreases from zero the neutral disturbances have two 
critical layers up to the point where the basic velocity is no longer non-monotonic. 
Then the two critical layers must merge, cease to remain asymptotically distinct and 
a fresh analysis is required. Since the principal thrust of the present work is to study 
resonant triads we do not explore this aspect any further here. However, we remark 
that these considerations show that neutrally stable linearised disturbances are pos- 
sible at propagation angles > tan -1 (^z*). There is therefore, a lower bound for this 
angle but no formal upper bound since from figure (3) we can see that neutral modes 
exist for /z. < /z < oo. This contrasts with the results of (12) where it was found that 
for linear disturbances the possible waveangles associated with upper branch neutral 
modes were bounded both above and below. 


Given the eigenvalues of the problem specified by (2.7), (2.8) fc (2.12) we sought 
resonant triad triples (oqj > Poj > fl()j ) of waves which are each neutrally stable at leading 
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order and which satisfy requirements (3.1). We found two such sets, in particular 

«01 =0.8514, /?oi = 0.264 8r n , fi 0 l = -3.137 x 10“ 2 , 

ao 2 =0.7343, /?02 = 0.0731r n , Cl Q2 = 1.12 x 1CT 2 , (3.3a) 

003=0.1171, /?03 = 0.1917r n , ft 03 = -4.25 x 10" 2 , 

and 

ooi=0.5847, Att = 0.2764r n , ft 0 i = -6.29 x 10" 2 , 

00 2 =0.3605, /?02 = 0.2699 r„, ft 0 2 = -7.16 x 10~ 2 , (3.36) 

003 = 0.2242, /? 03 = 6.46 x 10 " 3 r„, fio 3 = 8.76 x 1(T 3 . 

We do not claim that these are the only possible triads but the important result 
here is that resonant triples with the sought neutral stability properties do indeed exist 
within a three-dimensional boundary layer. The type of analytical and numerical work 
presented here can be applied to other basic velocity profiles. Indeed very similar 
results as to the forms of the neutral stability curves in oo//?o and Qq/(3q spaces 
appear for the attachment line profile and, of course, given the existence of triads 
for a simple rotating disc flow we have no reason to doubt that a similar conclusion 
holds for other much more complicated and more physically realistic three-dimensional 
boundary layers. 

This verification of the existence of resonant triads leads us on to an integral 
part of the current work- that of considering the disturbance structure of the three 
members of the triad when they are permitted to interact and evolve on a suitably 
slow timescale. 


§4. An asym ptotic description of the interaction mechanisms in a resonant 
triad of crossflow vortices. 

We have illustrated that within the context of three-dimensional boundary layer 
flows resonant triads comprising of crossflow vortices can arise. We now investigate 
an asymptotic description of how the members of the triad may evolve in time and 
space once the interaction processes are accounted for. To ensure that the structures 
of the critical layers in the flow remain as simple as possible for this study we assume 


12 


that the amplitude of each vortex is 0 ( 6 ) where <5 < 1 is sufficiently small such that 
all the the critical layers are linear in character. Then, through the nonlinear terms 
in the Navier Stokes equations, it is found that each pair of modes combine to drive 
the third mode via terms of size 0 ( 6 2 ) and this in turn suggests that the interaction 
mechanisms force the vortices to evolve on slow length and time scales of size 0(6). 

If again we relate our working to the rotating disc problem as a concrete example, 
suppose that we wish to describe the interaction of three linear vortices with leading 
order wavenumbers a o ; - , Po, and of frequencies floj ( j = 1, 2, 3) such that the resonance 
conditions (3.1) hold. Then we define 

Ej = exp [4y (y a ; - (r, e)dr + 9{3j (e) — tfi, (e)) j , (4.1a) 

^ r _ 


where 


ctj = aoj + eaij + . . . , P, = Poj + efilj + ■ • • , fl,- - flq,- + efliy + • • • > (4-16) 

and E\ = E 2 E 3 . We allow the disturbance vortices to develop on slow lengthscales 
ri = & ( r _ r n ), 6i = 66 and timescale ti = 6t , so that in the governing Navier Stokes 
equations we formally use the methods of multiple scales to replace the r, 6 and t 
derivatives according to the transformations 


d d _d_ d_ d ■ c g 

dri d 6 ~^dd d 6 i 

l + e ±. 

dt dt dt\ 


(4.2) 


As previously, we choose to examine vortices in the neighbourhood of some point 
r = r n and then, in the main inviscid region (I) where z = e 3 f , we are led to consider 
perturbation quantities of the forms 

3 


U = 6 


E (U 00j + ...)A j E j +c.c. 
L;=l 


+ 6 4 


E UiOjEj + ... 

Li =1 


[ 


1 . ,J. 


V = 6 j E (Voq# + ■ • •) A, Ej + c.c. j + 6 2 ^E v lQi E i+-- 


W = 6 


L- 

3 

E (W 00 j+...)AjEj+c.c. 


+ 


+ 


Li=l 


+ 6 4 


E WiOjEj + ... 


Li=i 


+ 


(4.3a) 


(4.36) 


(4.3c) 
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P = 6 


(-Poo, + • • •) AjEj + c.c. 

Lj=l 


+ 6 : 


S PlO; Ej + 


Li =1 


+ .... 


(4.3d) 


In (4.3) the scaled amplitudes A, are functions of the slow scales r\, 9\ and fi 
and U 00j , V 00 ; , • ■ • > PlO* > • • • are functions of £ alone. Substituting (4.2) and (4.3) into 
the Navier Stokes equations and comparing coefficients of Ej we find, as in section 2, 
that at leading orders Wqqj satisfies the Rayleigh equation 


v* (w 00h -70*^00*) - &0k W 00k =0, k = 1,2,3, 


(4.4) 


where D jk = a jk r n n+p jk d~n jk and jf k = aJ k +(P? k /r 2 ). The functions U 00} and V 0 0; 
may be found in terms of l+oo ; - by suitable manipulations of the governing equations, 
which are precisely (2.6) with a 0 ,/?o,^0 replaced by a 0j , (3 0j , and Uq,V 0 ,W 0 ,Pq 

by Uooj , Vooj , Wooj , Poo/ ■ 

^From here onward we concentrate on deriving the evolution equation for the 
vortex proportional to E\. Given the obvious symmetry in the problem, together 
with the generalised way in which the amplitude equations will be derived, once the 
governing equation has been obtained for one mode it is formally straightforward 
to obtain the corresponding equations for the other two modes. At 0(6 2 ) in the 
Navier Stokes equations we find that the function l+ioi in (4.3c) satisfies the non- 
homogeneous Rayleigh equation 

ffoi ( 1+101 — 701 l+10l) - ffoi^lOl = 

_ iWm - 2,^01^001 («01^ + ^^) + I.A 2 A 3 , 

(4.5a) 

where 


I* — * [(/?03<ar01 + a 02^0l) 002 + (/?02 Q! 01 + O(03A)l) F 002 ^ 003 + 0(01^002^003+ 

v r n r„ 

2 * i 

^r+002+003| + 1+003 ( 001^002 + — +002^ + 1+002 ( 001^003 + — +003) 

r n J L ' r n / V r n / 

— *701 [* ( 002^003 1+002 + 0(031+003^002 + — +0031+002 + — +OO2I+OO3) 

^ ? n fix * 


+I+003I+002 + I+002I+003 


(4.56) 
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and 


(4.5c) 


( d d d \ 

W r " B ftT + W' 

The appropriate boundary conditions on the function Woi very similar to 
those for Wool in (4.4). In particular, to avoid unbounded growth as £ — ► oo requires 
that |Wioi(oo)| < oo and for no disturbance at the surface of the disc, Wioi(0+) = 
0. Additionally, it is clear that appropriate jump conditions need to be imposed at 
the position of the critical layer (s). For ease of notation, suppose that the mode 
proportional to Ej has associated critical layers located at and note that there 
could well be more than one critical point for each mode. Inspection of (4.5) suggests 
that singularities in the solution Wioi(0 may be expected at all the critical layers 
for j = 1,2,3, not just at the positions £ l . 

The form of the solution of (4.5a) in the vicinity of the critical layers at may 
be one of two distinct types. Close to the position £i we can easily show that the 
right-hand side of (4.5a) is asymptotically proportional to l/(£ — £i) and from this we 
can deduce that for 0 < £ - £i < 1, 

W m = Xo log(£ - &) + 0(1) + . . . , (4.6a) 


where Xo is a constant which may be found analytically. Close to the other critical 
layers, £,• (;' = 2,3), the behaviour of Wioi is more complicated. Then we obtain, from 
(2.6) k (2.8), that 




* 0 ; , Xu ^ X 2j 

(€-$)* ( e -£) 2 tt - 5 ) 


+ 


(4.66) 


and then 

Wm ~ TT^r\ + Dl * - Si ) + D 2j + 0 3i (£ - ) ln(e -£,) + ..., (4.6c) 

where the constants X kj , D kj are given in the Appendix. The solutions (4.6) are valid 
for £ > £; and clearly need to be modified for the region £ < This requires careful 
scrutiny of the critical layer properties in the vicinity of f in much the same way as 

was outlined in §2. The main difference now is that the algebraic manipulation needed 
within the critical layer region becomes much more complex due to the presence of 
all three modes. Rather than give details of the inevitably lengthy (but conceptually 
fairly straightforward) critical layer analysis, we elect to merely state that the result 
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of this work is as is to be expected from the classical theory: that is that we obtain a 
‘?V phase shift in the logarithmic terms as characterised by (2.12). This then enables 
us to write down the behaviour of W m for 0 < £• - £ < 1 and thence to continue the 
solution ITioi into the region £ < . 

To derive the evolution equation for the amplitude A j we follow the usual process 
of multiplying (4.5a) by PTooi and integrating over the range [0, oo) , paying particular 
attention to the behaviour of the integrals close to the critical layers. Recalling the 
solutions (4.6) and the phase jumps across the critical layers, combining the integration 
with the boundary conditions imposed at £ = 0 and at £ = leads to the desired 
amplitude equation. After simplification, and substituting the definition of C. given by 
(4.5c), we obtain 


( 27tBqi 2 _ „- TTy („/ \ V T n U 2 , 

finite part |i?oil| 001 




nr. 


• n u f Bqi2 Bqiz\ 2 *^ 2r n 0^012^001 


~Sf 


C (w W ™ 


T n uV, 


n 'fr2 01 ^001 j + 2 i J aoir n WQ 01 dt 

u 0 l * ^O 


Oil 

OO 


^001 + 

). 


01 


|-Soil|5on 


dA i 
J dri 


+ 




7TV 


\ W 001 + 


7rBoi2VW& 


001 


7TC 


finite part |5 0 n| " Wi ’ |J?oil|£oil 
( B oi2 5 0 i 3 \ 2 

ym: b^> w °° i mr. — J - 


l^on | x -Son jBoh 

+l h ) d * + 2i L ^ lW ooi d ^ 


do 1 


(4.7) 



n 

^012 

■Sols'! 

2 l7r b $ 12 W'ooi^ 

[V 

|S 0 n| 

' -Soil 

■Soil'' 

^001 ^ J 

■Don 7 


_j / !i|iv2 

K Wl 


U 01 w 2 At 

m w md( \ i>h 


+ 


V* \ ilE 2 il m n n « *1 [°° LW 001 

^ 1 | Bqj 1 1 {Djl9jl D *9j0)\ -A - 

u 


i= 2 




d{ 


A 2 A 3 = 0. 
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In this amplitude equation we have defined 


B 0jlt = (r n ao*G (k) (0 + Po *® ( °(0) . ( 4 - 8 ) 


and have supposed that 

00 

*001 = E (*-£)*. 

i=0 

close to the critical layers £j- ( j = 2, 3) at which this function does not have singularities. 
Otherwise, the notation is as previously defined. 

In (4.7) the expression / denotes the Hadamard (finite) part of the corresponding 
integral. Some of the integrals which are contained in the coefficients of dA\/dr \ , 
dAi/ddi and dAi/dti are singular at fi: the integral in the coefficient of A 2 A 3 is 
singular at the other critical points as well. 

We notice that if the various integrals and summations were to be evaluated nu- 
merically and the formulae (4.7) rearranged appropriately, we could obtain an evolution 
equation for the E\ mode which assumes the form; 

+ A13A2A3 = 0, (4.9a) 

dt\ dr 1 do\ 

where the coefficients An,Ai2>Ai3 arc some constants. By repeating the analysis out- 
lined above we can, in an identical way, derive evolution equations for the other two 
modes of the types 

^ + A 2 i^2 + A 22 ^ + = 0, (4.96) 

dt 1 dr 1 o&i 

+ Ml—- + A 3 2^p + AsaAiAJ = 0, (4.9c) 

dt\ dr 1 do 1 

where an asterisk on a quantity denotes the complex conjugate of that quantity. 

Clearly, the numerical work required to actually evaluate the coefficients in (4.9) 
for any particular basic flow is nontrivial. However, to illustrate the type of evolution 
system which may be obtained we considered the computational procedure relating to 
the resonant triple given by (3.3a). For each member of the triad we computed the 
eigenfunction Woo, by solving the Rayleigh equation (2.7 a) subject to boundary con- 
ditions (2.76, c) and jump conditions (2.12). The corresponding radial and azimuthal 
velocity fields IToo,- , Voo, may then be determined in terms of the functions Uoj and 
Woo,- by eliminating the pressure from (2.66, c). Since the problem (2.7) is linear we 


17 


needed to fix some normalising constraint and this was chosen such that as £ ► oo, 

Woo, was a real, positive function and that max{0 < f < oo : 1^00*1 = l}. 

At this stage the coefficients A ; i and A, 2 in ( 4 . 9 ) were evaluated. The integrals 
involved in finding these unknowns are singular only at £• and the analytic nature 
of these singularities may be deduced in a routine manner. To find the finite parts 
of the integrals across the positions , we subtracted the singular behaviour of the 
integrand in the neighbourhood of £ and integrated the remaining (necessarily well 
behaved) function across fj. By utilising the formula given by Lighthill in ( 26 ) it 
was possible to deduce the Hadamard part of the subtracted singularity function and 
thence determine the value of the finite part of the original integral. The method 
worked well in practice and the corresponding results are given below. 

We turn now to address the much more complicated task of determining the 
interaction coefficients Aj,- 3 and again, as a concrete example, concentrate on the de- 
termination of A13 for the method for calculating the other interaction coefficients 
follows in a very similar way. Then the function I, is as defined in ( 4 . 56 ) and was 
computed by using the information previously derived on Uoqj , V 00j and W 00j and 
detailed above. It can be shown by straightforward analysis and using (4.66) that the 
integrand (/, Wooi/C^oi) has singularities at each critical layer and to calculate this 
integral a technique was employed which is very similar to that already described. In 
the vicinity of each critical layer the singular part of the integrand was removed as 
before, except now we found that the remaining function was very erratic near the 
critical point. This is due to the fact that the integrand can only be determined nu- 
merically at the outset and the process involves the estimation of several derivatives of 
Uooj , Voo; which is a probable source of numerical error. On making the usual checks 
for numerical stability, the determination of the finite part of the integral was found 
to be strongly influenced by the choice of the intervals around the critical points in 
which the singular behaviour was removed. This is obviously unreasonable and a more 
satisfactory procedure was found in which the integrand was written as the sum of 
the singular part and an interpolated polynomial in the vicinity of the critical points. 
Even so, in this case there was found to be some variation in the determined value 
of the coefficient A 13, but it has to be remembered that this coefficient depends on a 
calculation in which four separate critical layers are traversed with a strong singularity 
needing to be accounted for at each critical layer. Consequently, the difficulties en- 
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countered in obtaining accurate values for the interaction coefficients in the evolution 
equations (4.9) are not surprising. 

In this paper we have managed to obtain approximate values for the interaction 
coefficients. What we have shown is that determination of the interaction coefficients 
for the flow structure considered here is possible in practice, and, with the use of more 
sophisticated numerical methods, there is no reason to doubt that the accuracy of the 
coefficients could be markedly improved. 

We now describe some properties of the solutions of the evolution equations and 
conclude with some discussion. 
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§5. Results and Discussion. 


The numerical work described above leads to evolution equations which assume 
the forms 

^ = r n (-0.114, 4.09 X 10“ 3 )^ + (0.305,0.111)^ + (2.0, 0.6 )A 2 A 3 , (5.1a) 

ot i or i Ou\ 

^ = r„ (7.52 X 10~ 2 , —7.74 x 10" 2 )^ + (0.265,0.117)^ + (0.15, -0.15)AiA^, 
ot i or i ov i 

(5.16) 

^ = r„(— 2.4 x l0- 2 ,-0.14)^ + (0.14, 0.13)|^ + (0.3, 0.4)AiA£, (5.1c) 

for the resonant triple given by (3.3a) and (4.3). As already mentioned, the effect 
of having four critical layers within the flow makes the calculation of the interaction 
coefficients very sensitive and we believe that the values of these quantities quoted in 
(5.1) are correct to within about ten per cent. 

Naturally, the evolution of the wave amplitudes Ai, A 2 and A 3 is very dependent 
upon the initial conditions and the behaviours of these quantities may, in general, only 
be determined by a numerical approach. We do not study this facet of the problem here 
but, instead, draw an analogy between the forms of (5.1) and those derived by Smith & 
Stewart (22) in their investigation of resonant triad nonlinear interaction in boundary 
layer transition. The evolution equations derived by these authors were a simpler form 
of ( 5 . 1 ), (effectively with djdr\ and d/d9\ replaced by constants) with the interaction 
coefficients purely imaginary quantities. Smith & Stewart integrated their equations 
using a predictor— corrector formulation to march forward in time and the resulting 
evolution of the wave amplitudes took a very complicated form. These amplitudes 
generally increased with time although this trend was modified by the superposition 
of a very short period, large oscillation which had the overall effect of giving the 
function | Aj (ti) | a very spiky appearance. Similar types of evolution equations have 
been derived by Avis (27), who was concerned with resonant triads in axisymmetric 
flow on a cylindrical geometry, and by Bowles (28), who investigated triads within 
transonic flow. These authors obtained numerical solutions for the development of the 
mode amplitudes, and both found behaviour qualitative to that described in ( 22 ). 

Very little is known about particular solutions of (5.1), although some work on 
related equations has been studied by Craik ( 1 , 29, 30). No general analytic solutions 
for these equations have been found. However, an important property of the system 
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(5.1) is under what conditions these equations admit singularities in the solutions A j, 
0’ = 1,2,3), either in a finite or an infinite time. Wang ( 31 ) has derived bounds for 
the solutions of three-wave interaction systems with arbitrary coupling coefficients. 
He showed that if 


da o 
dt 


= iwQQQ + n oaia2, 


da\ , da 2 . , 

— = iwiai + Viaoa 2 , — = iu 2 a 2 + H 2 aoa\, (5.2) 


with initial data a<(0) = a<o, * = 1,2,3: then explosive instabilities (instabilities that 
become unbounded eventually, possibly after infinite time) cannot exist if 


a < 0 and | ja(0)|| < V\a\/0, (5.3a) 

where 

a = min{ max{Aoo, A 12 , A21}, max{Aoi, A10, A22}, 

max{Aoo, A01, A22}, max{Aoi, A12, A21}, ( 5 . 36 ) 

max{Ao2, A10, A21}, max{Ao2, An, A20}}, 

P = \ max{|^ |}, and A \ - 21m (w*) . (5.3c) 

Here ||a|| = || (a 0 ,ai,a 2 ) || = (E.^oaiarJ . 

If we allow the disturbances proportional to (4.1) to be neutrally stable in space 
co-ordinates r\ and d\ we can substitute 

ia 1 

and ifi 1 , (5.4) 

for d/dr 1 and d/d9\ in (5.1) with c*i and 0\ real. Then using the results (5.3) we 
can find the region in ( a\,P \ ) space for which explosive instabilities of the system 
(5.1) are impossible, at least for sufficiently small initial disturbance amplitudes. This 
information is illustrated in Figure (4) where we indicate the minimum size of the initial 
amplitude required in order to ensure that the possibility of explosive breakdown is 
maintained. Broadly speaking, we can conclude from figure (4) that once we are 
located in that part of (ai,/?i) space at which a finite initial amplitude is required to 
enable explosive breakdown then by decreasing the radial wavenumber ai and/or by 
increasing the azimuthal wavenumber there needs to be a corresponding rise in the 
initial disturbance amplitude if the breakdown mechanism is to remain possible. 
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It should be noted that the point corresponding to the choice ai = = 0 

lies in the zone for which the possibility of explosive instability cannot be discounted 
irrespective of the initial conditions. Then the governing equations reduce to the forms 


-j~ - = (2.0, 0,6)A2j4.3, — — = (0.15, — 0.15)AiA 3, -jr^ = (0.3, 0.4 )Ai^2- 

at i ati at i 

(5.5) 

These are, of course, the evolution equations obtained by allowing the disturbances 
to evolve temporally whilst suppressing spatial dependence. It is a straightforward task 
to write i>i = 62 = A 2 , 63 = A3 and to normalise so that the system (5.5) may be 


rewritten as 


-TT = Sib*b* k , (i,j,k= 1,2,3 cyclically), 
aci 


where the interaction coefficients = exp (ifa). Craik (29, 30) has shown that ‘phase- 
locked’ solutions of systems analogous to (5.6) exist where © = Q\ + 62 + #3 is constant 
and &i = phase(bi). Then 0 must be a root of 


Y tan(<^ — 0) = 0. 


In this case the existence of roots of this equation is guaranteed and then, for ap- 
propriate initial conditions, the solutions can exhibit a finite time singularity. It 
has to be remembered that in practice once the disturbance amplitudes become too 
large however, then the assumptions concerning the smallness of the amplitude within 
definitions (4.3) cease to remain valid and a revised asymptotic analysis is required. 

Apart from the results relating to explosive instabilities described above, very little 
is known about solutions of system (5.2) for general parameter values. In the course 
of an investigation into interaction processes between crossflow vortices in practically 
important cases one would have to resort to extensive numerical calculations to provide 
a f ull description of the characteristics of the solution of the evolution equations. This 
provides scope for extending the elementary ideas presented here to the potentially 
more important interactions in experimentally generated boundary layers. 

The derivation of the evolution equations relating to the second resonant triad as 
given by (3.3 b) was also considered. However, for this triad, the mode proportional to 
E 3 has a critical layer very close to the wall. In practice, it was found that the numerical 
scheme described in §4 was too crude to accurately model the singular behaviour close 
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to the wall to provide satisfactory convergence. This discovery points to the need to 
develop improved numerical procedures to compute the interaction coefficients. 

The asymptotic arguments behind the derivation of the amplitude equations for 
modes with spatial and temporal dependences as in (4.1), (4.2) relied on the assump- 
tion that the amplitudes of the disturbances were 0{6) across the majority of the 
boundary layer, where <5 < 1. This condition was introduced primarily to ensure that 
the viscous critical layers within the flow assumed their simplest possible forms: i.e. 
linear critical layers encompassing the usual ‘tV phase shift. We can formally relax 
this requirement by making the following observations. Recall that 

Ej sexp \\(J' a^r^dr + ePiW -&,('))] . (* la ) 

and suppose that, according to linear theory, the mode proportional to Ej is neutrally 

stable when k _ 

aj =a Oj + totij + + € a hj + • ■ ■ > 

/ 3 ;■ =(3oj + efilj + + t k Pkj + •■■■> ( 5 - 7 ) 

Ctj =Ooj + eftij + + + — 

Suppose we allow the disturbance to evolve on slow space and time scales given by 
(4.2) when we replace 

6 by e* (^- + *«*,)» ( 5 - 8 ) 

dr i 'dr i / 

and make similar replacements for 6{df89i) and 6(0/0^). Then the analysis of the 
interaction mechanism will proceed in a very similar way as before, so long as the 
critical layers remain linear in form. Following the work of ( 8 , 20 ), we know that 
the critical layers appearing in this type of disturbance structure remain Unear in 
character until the disturbance amplitude reaches the size 0(e 2 ) at which point the 
critical layer configuration is modified due to the presence of (relatively) large mean 
flow corrections and the calculation of the phase shift across the layer becomes a 
numerical task. However, we can see that the derivation of the evolution triad of 
equations (4.7) may be trivially modified to account for the interaction between three 
modes of finite amplitudes, so long as these ampUtudes remain less than 0(e 2 ). 

To conclude, we beUeve that the major finding of our work here is that triads 
of Unearly neutrally stable crossflow vortices can exist in three-dimensional boundary 
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layers; further these triads are such that we have a resonant set of modes at lead- 
ing orders in wavenumber/ wavespeed expansions based on the asymptotically large 
Reynolds number. Further, we have developed a fully rational description of the in- 
teraction processes which promises to reveal the degree to which these mechanisms 
can alter the stability characteristics of the flow. We have concentrated on a special 
three-dimensional flow as an illustrative example but the results obtained here suggest 
that the interaction process will be operative in more practically oriented boundary 
layer flows. Our work can be regarded as a first stage in a very complicated field: 
obvious generalisations could be tackled to explain situations in which, for example, 
we have critical layers of other than the simple linear kind; or in which the modes 
have slight detuning. In this context we note that Avis (27) has made a careful study 
of the effects of detuning on resonant triads in his axisymmetric flows. The analysis 
adopted here allows nonparallel effects to be accounted for without formal difficulties- 
these effects can play significant roles in practical applications. The largest drawback 
of our approach is that the computation of the interaction coefficients appears to be 
not straightforward, especially when the modes have critical points close together or 
when such critical points are near the bounding surface. We feel that to be able to use 
our approach to its best advantage further study needs to be devoted to developing 
more powerful and sophisticated techniques for resolution of the numerical problem 
as opposed to the fairly rudimentary ideas employed here. Overall however, we have 
here the first asymptotic description by which the importance of interaction mecha- 
nisms upon a resonant triad of linear neutrally stable crossflow vortices may be truly 
assessed. 
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Appendix. 


Here we briefly describe the derivation of the constants {A*,' , D kj } appearing in 
(4.6). We shall, for convenience, determine these quantities when j = 3 for then we 
may easily modify the results for j = 2. 

With I, as defined in (4.56), we suppose that close to the critical layer £3 the 
eigenfunction W 003 develops according to 

W 003 ~ dl 3 + 0 ( (£ - 6) In K - 61) ■ ( A1 ) 

Further, in this region the mode proportional to E 2 has a regular Taylor expansion 
and so 

(^002,V002,W002) ~ (Qu02) Qv02i Qw 02) + U ~ 6) (Qtt02) Qv02> Qw 02 ) + •••» (^ 2 ) 


where the quantities Q|*q 2 } are constants. Finally, we know that the behaviours of 
Uqq$ and Vqq 3 can be obtained by suitable manipulations of the basic equations (2.6). 
Substituting this information into (4.56) yields the behaviour 


A 03 -^13 . -^23 

(* - 6) 3 (£ - 6) 2 (€ - 6) 


(A3) 


where 


X03 = _ 2 *^™ 02 ^ 031 ^ (c*0lA)3 - «03/?0l) . 
7TP03-D031 


Xu = — (ao 3 / 5 oi - ctaipaa) [(/?oi + £03) Q v 02 + (<*01 + <*03) < 3 * 02 ] 

Trp 03-0031 


ig| 3 Ho32^01lQ u ,02 ^|-0Q3ll 


(A4) 


(A5) 


jdT 

-D031 


7r/?03-Bo31 


(«0l/?03 - «03/?0l) Qt»02> 


and 


= _ ig,o27§i|Bo3ij {ao2 ^ _ a03/)02) j_ 

5g 31 7T/?03-D031 

g °pf 032 ~ [( a 0 lQ u 02 + /? 0 lQ„ 02 ) -Boil + (a: 03^01 - « 0 l/? 03 ) {Q%02% 
B 03l 


Qv02r n ®i 3 )] ■ 


(A6) 


Here 


7rao3-Bo315(L _ ^®e^£03 
703l-B03l| 1-0031 1 
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Then we can show that 


WiOl 


D 03 


_ ~Q + -°13 In (f ~ fa) + T>32 + I>33 (f - f 3 ) In (f - f 3 ) + 


(A7) 


for £ > £3 with 


&03 


*03 

2i?oio ’ 


Da = 2 -^d 0 3 - * 13 


-B 010 


•Boio’ 


and 


•°33 = t^I>13 + 701-^03 + "^ 33 


Brno "” ' ' U1 “ U,, 1 £010 

The quantities £> 13 , D 33 appear in the evolution equation (4.7) for the mode pro- 
portional to E\. Of course, we need to obtain quantities equivalent to those given 
above in order to traverse the critical layer(s) £ 2 , but these can be deduced merely by 
interchanging the subscripts ‘2’ and ‘3’ in the relevant formulae above. 
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Figure (1). 

Flow structure for the Rayleigh modes described in §2. The inviscid layer I is of 
thickness 0(e 3 ) (the same depth as the boundary layer and of the same size as the 
horizontal scale for the structure). In addition, we have a viscous wall layer {II) of 
depth 0(e 9/2 ) and a critical layer (III) of thickness 0(e 4 ) at those positions where 
the ‘effective’ basic velocity profile Cq (see (2.7)) vanishes. 










Figure (3). 

Solutions of the eigenproblem (2.7) for various values of the waveangle parameter A*- 
a) Radial wavenumber ao as a function of the azimuthal wavenumber /Jo- 
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is sufficiently large: the minimum values of A(0) needed for this are indicated on the 
solid contours. 
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